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We study the geometric phase of the ground state in the extended quantum compass model in 
presence of a transverse field. The exact solution is obtained by using the Jordan- Wigner transfor- 
mation which maps the Hamiltonian on a fermionic system and applying the Fourier transformation 
to realize the diagonalization and an analytical expression for the ground state and geometric phase 
in the momentum space. Furthermore, the scaling behavior of the extremum of geometric phase 
and its universality are investigated due to the divergence of geometric phase at the critical point. 
We also study the fidelity of the ground state to confirm the results obtained using the geometric 
phase. 
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I. INTRODUCTION 



Quantum phase transitions (QPTs) are essential phe- 
nomena in quantum many-body correlated systems 
which is induced by the ground-state (GS) transition. 
This can be driven by tuning an external parameters at 
zero temperature, for instance, the magnetic field or the 
amount of disorder. Traditionally, QPTs are described 
in terms of order parameter and symmetry breaking by 
the Landau-Ginzburg paradigm, the problem in these 
traditional ways is that there is no general method to 
find the order parameter for a common system. In the 
last few years a big effort has been assigned to the anal- 
ysis of the QPTs from the perspective of Quantum In- 
formation [T]-^. Entanglement and fidelity have been 
accepted as new notion to characterize quantum phase 
transitions. Entanglement, referring to quantum correla- 
tions between subsystems is a good indicator of quantum 
phase transitions, because the correlation length diverge 
at the quantum critical points [IHl]- The fidelity which 
is a measure of distance between quantum states, could 
be also a nice tool to study the drastic change in the 
ground states in quantum phase transitions [TD] . In view 
of some difficulties [TT], attention has been transferred 
to include other, potentially related, means of charac- 
terizing QPTs lOj. One such approach centers around 
the concept of geometric phase (GP). Geometric phases 
which display the curvature of the Hilbert space, and 
especially has direct relation to the property of the de- 
generacy in quantum systems, were introduced in quan- 
tum mechanics by Berry in 1984 [T^]. Thus the GP is 
a powerful tool to find out the quantum phase transi- 
tions [13]. The attention has also increased due to their 
applicability in quantum-information processing |14j . In 
other words, GP becomes additive for product states of 
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composite systems since the uncorrelated subsystems set 
independent geometric phase factors. However, due to 
their global properties, the GP is propitious to make 
fault tolerant quantum gates. In this point of view, many 
physical systems have been studied to realize geometric 
quantum computation, such as NMR (Nuclear Magnetic 
Resonance), Josephson junction, ion trap and semicon- 
ductor quantum dots. The quantum computation scheme 
for the geometric phase has been proposed based on the 
Abelian or non-Abelian geometric phase, in which geo- 
metric phase has been used to be robust against faults in 
the presence of some kind of external noise due to the ge- 
ometric nature of Berry phase. It was therefore seen that 
interactions play an important role for the understanding 
of some specific operations. 

For the first time A. CaroUo at ai, show the GP con- 
nection to the quantum phase transitions, where they 
have been shown that the geometric phase of the spin 
chains diverge at the critical point [15]. The critical ex- 
ponents were evaluated from the scaling behavior of ge- 
ometric phases has been obtained by S. L. Zhu at al, 
[TB]. Then the geometric phase have been considered as 
a topological test to manifest quantum phase transitions 

In the present work, we shall use the geometric phase 
and fidelity of the ground state to detect the QPTs for 
the one-dimensional extended quantum compass model 
pni [m [55] in a transverse filed to confirm the results ob- 
tained in our previous works [THl [TS] . This model cover 
a group of well-known spin models as its special cases 
and shows a rich phase diagram. It is worthy to men- 
tion that not only this inhomogeneous system exhibits 
rich phase diagram, it would be interesting to investi- 
gate whether the geometric phase is able to determine 
the quantum phase transition in these more complicated 
systems. However, an explicit analytical expression of 
the ground state, which is necessary for studying the ge- 
ometric phase and fidelity, has not been calculated so far. 
The exact solution is obtained by using Jordan- Wigner 
(JW) transformation and because of nice scaling proper- 
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ties of GP, phase transition can be pinpointed from small 
systems with considerable accuracy. 



Hamiltonian Eq. ^ can be written in the diagonal 
block form 



II. HAMILTONIAN AND EXACT SOLUTION 
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where Ji and J2 are the odd bonds exchange couplings, 
Li is the even bond exchange coupling and N = 2N' 
is the number of spins. We assume periodic bound- 
ary conditions (ground state for periodic and antiperi- 
odic boundary conditions are identical in the thermo- 
dynamic limit and the essential features in finite size 
are also not altered qualitatively). This model em- 
braces a group of other familiar spin models as its spe- 
cial cases, such as quantum Ising model in a transverse 
field for J2 = L2 = 0, the transverse field XY model for 
Ji = Li, J2 = L2 and the transverse field XX model for 
Ji — J2 = Li = L2. The above Hamiltonian (Eq. Q) 
can be exactly diagonalized by standard Jordan- Wigner 
transformation ilSi ilOi |22i as defined below, 
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which transforms spins into fermion operators Cj. 

The crucial step is to define independent Majorana 
fermions [231 El] at site n, = C2n-i and = C2n- 
This can be regarded as quasiparticles' spin or as split- 
ting the chain into bi-atomic elementary cells |23| . 

Substituting for crj, cr| and cr| (j = 2n, 2n— 1) in terms 
of Majorana fermions with antiperiodic boundary condi- 
tion (subspace with even number of fermions) followed 
by a Fourier transformation, Hamiltonian Eq. ([I]) (apart 
from additive constant), can be written as 



k 



where J = (Ji - J2) - (Li - i2)e*^ i = (Ji + J2) + 
(Li + L2)e'^ and k = (j = 1, 3, • • • , A^' - 1). 

By grouping together terms with k and —k, the Hamil- 
tonian is decoupled into a sum of independent terms act- 
ing in the 4-dimensional Hilbert spaces generated by k 
and —k (7J+ = 0j,>o ^k)^ ^he other word [Hk,Hk'] ~ 
in which 



^^k Cfe - Cj,C_j, - Lj C^C-f, 
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where ^ = {4^c%4,c?X) and 
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By using the element of the new vector F'J 



(7fci 7_fe, 7fc: 7?fe) could be described by unitary trans- 
formation Vj. = Uk^k (see Appendix A) the matrix A{k) 
can be diagonalized easily and we find the Hamiltonian 
Eq. ([2]) in a diagonal form. 
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where E? — \/ a + \/h and E\ ~ \J a — Vb, in which 



a 
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Ah^ + \J\^ + \L\'', 
(16/i2 + 2|J|2)|i|2 



The ground state (Eg) and the lowest excited state 
{Ee) energies are obtained from 



Er 



It is straightforward to show that the energy gap 
vanishes at ho = ^ (Ji + L2){J2 + Li) and h^^ = 
\J {Ji + L2){J2 — Li) in the thermodynamic limit. 

So, the quantum phase transition (QPT) which could 
be driven by the transverse- field, depending on exchange 
couplings, occurs at /iq and /i^r- 

By tedious calculation on the unitary transformation 
the unnormalized ground state has been written in the 
vacuum /cth mode of 4 ^'^d c^, 

\G{h)) = n ['^ilO) + "2 cf cl\.|0) + z;3 cf c^JjO) 



+^6 clVf C^JJO)], 



(5) 



where Vi, (z = 1, • • • , 6) is functions of the coupling con- 
stant which are given in Appendix B. We should point 
out that there is a flaw in Ref. [25] in the basic diago- 
nalization results of the Hamiltonian which is the same 
as the present model in k-space. 
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FIG. 1: (Color online) Two-dimensional phase diagram of the 
extended quantum compass model in a transverse magnetic 
field for zero temperature (Ji = 1). 
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FIG. 2: (Color online) Zero-temperature phase diagram of the 
extended quantum compass model in the transverse magnetic 
field for Jl = -1. 



III. PHASE DIAGRAM 

It is useful to recall the zero-temperature phase di- 
agram of the extended compass model in homogenous 
transverse field, see Figs, ([l] [2]) [HI [19] (For simplicity 
we take Li = 1, L2 = 0). In our recent work, we have 
investigated the phase diagram of the extended quantum 
compass model in homogenous transverse field by use of 
the gap analysis and universality of derivative of the cor- 
relation functions [THj. However the phase diagram [25] . 
quantum correlation [27] , bipartite entanglement [55] and 
fidelity [29] of this model has been studied numerically 
using the exact diagonalization method and infinite time- 
evolving block decimation. This model is always gapful 
except at the critical surfaces where the energy gap dis- 
appears. We have obtained the analytic expressions for 
all critical fields which drive quantum phase transitions 
(QPT) as a function of exchange couplings 15], 

ho = y/MJ2 + Li), K = \/Ji{J2-Li) for L2 = 0. 
There are four gapped phases in the exchange couplings' 
space: 

• Region (I) Ji = 1, < J2 < 1: In this region for 
small magnetic field {h < ho) the ground state is 
in the spin-flop phase (Fig. ([!])). 

• Region (II) Ji — 1, J2 > 1: In this case there is 
antiparallel ordering of spin y component on odd 
bonds for h < h^. In this region tuning the mag- 
netic field forces the system goes into a spin-flop 
phase (region (I)) for hj^ < h < ho- ior h > ho the 
system fall into region (III) (Fig. ([I])). 

• Region (III) Ji = ±1: In this region the ground 
state is the ferromagnetically polarized state along 
the magnetic field (Figs. ([l]),(|2|). 



• Region (IV) Ji = — 1, J2 < 1: In this region the 
ground state is in the strip antiferromagnetic (SAF) 
phase (Fig. 

The fidelity and geometric phase of the ground state 
properties could reflect these different zero-temperature 
regions, in particular the quantum phase transitions at 
ho and h^^. 

IV. GEOMETRIC PHASE AND FIDELITY 

In order to investigate the geometric phase in this sys- 
tem, we use a new family of Hamiltonians that can be 
described by applying a rotation of ip around the z direc- 
tion to each spin as in the method in Refs. [15] and [16], 
i.e., 

N 

H4, = g^Hg^, ff0 = exp(i(/)cr|/2). 

The energy spectrum does not depend on the angle cj), 
and the critical behavior is independent of </>. The cor- 
responding ground state of the new Hamiltonian may be 
obtained as \Gc/,{h)) — g^\G{h)) which is described by 

= n t^ilO) + ^2e2'^cf cl\|0) + e^'^acf c^IlO) 

k 

The geometric phase of the ground state will be accu- 
mulated when the system finish a cyclic evolution, corre- 
sponding to varying the angle (j) h'om to tt, 

^^-^ j\G,{h)\^\G^{h)) 



FIG. 3: (Color online) The geometric phase of quantum com- 
peiss model in a transverse field as a functions of transverse 
field and J2 for Ji = 1 and iV = fOO. 



FIG. 4: (Color online) The geometric phase of the model 
versus h and J2 for Ji = — 1 and A*' — 100. 



and is found to be 



/3s = 
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By this expression, we could study the geometric phase 
of the ground state and its derivation to investigate the 
QPTs of the model. 

However a sudden drop of the fidehty, caused by the 
ground-state level-crossing is too obvious to be interest- 
ing enough, to detect the existence of the QPTs. Those 
ground-state wavefunctions are interested in which are 
differentiable in parameter space. Therefore, the overlap 
between two ground states at h and h + 6h can be defined 



F{h,Sh) ^ \{G{h)\G{h + 6h))\, 



(7) 



in which Sh is a small deviation in a transverse field. 

Three-dimensional panorama of the geometric phase of 
the model has been plotted as a functions of a transverse 
field and J2 in Figs. (|3| and ^ for Ji = 1 and Ji = -1 
respectively, and system size set as = 100. In Fig. 
([3]) it has been observed that the geometric phase has 
its minimum value for < /i < /i^ and J2 > 1, which 
corresponds to the region (II) in the phase diagram (Fig. 
([1])). It manifests that the geometric phase stay quite 
unchanged in region (II) . The geometric phase undergo a 
strong qualitative change in spin- flop phase (region (I)) 
and enhances by increasing the magnetic field. In region 
(III) which corresponds to the saturate ferromagnetic 
phase {h > ho) increasing the magnetic field saturates 
the geometric phase too. The strip antiferromagnetic 
phase (region (IV) in Fig. (I2I) has been distinguished 



from the saturate ferromagnetic phase by the concave 
surface in Fig. ([4]). However, Figs. ^ and ([4| survives 
that the geometric phase enhances with increasing the 
magnetic field. We should point out that the two criti- 
cal lines {ho{J2, h) and hT^{J2, h)) at which the quantum 
phase transitions occur can be described by two assumed 
lines on the convex parts of the surfaces in Figs. ([3| and 
([4]). It is expected that the derivative of the geometric 
phase probes the QPTs. Then it will be helpful to study 
the first derivative of the GP with respect to the mag- 
netic field, three dimensional plot of the derivative of the 
GP with respect to h has been shown in Figs. ^ and 
([6]) as a functions of h and the Hamiltonian parameter 
J2 for N = 100. There are clear peaks for the derivative 
of geometric phase near the critical lines hQ{J2,h) and 
h-K{J2, h) for finite size lattice. Although there is no real 
divergence for finite lattice size, but the curves exhibit 
marked anomalies with height of peak increasing with 
the system size. However, similar to the geometric phase, 
the notation of fidelity which also extract from the field 
of quantum information theory can be used to character- 
ize the QPTs. The basic idea is that near a QPTs point 
there is a intensive enhancement in the degree of distin- 
guishability between two ground states, corresponding to 
different values of the parameter space which defines the 
Hamiltonian. This distinguishability can be specified by 
the fidelity, which for pure states reduces to the ampli- 
tude of inner product or overlap. The behavior can be 
ascribed to a dramatic change in the structure of the 
ground state of the system during the quantum phase 
transition. Therefore, it should contain all the informa- 
tion that describes QPTs and topological order. 

The ground sate fidelity of the model has been depicted 
in Figs. ([7| and (|8| for Ji = 1 and Ji = —1, where we 
have set iV = 100 and 6h = 0.001. Obviously, there is 
a sudden drop in the ground state fidelity at the QPTs 
lines. Comparing the divergence of the geometric phase 
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FIG. 5: (Color online) Three-dimensional panorama of the 
derivative dPg/dh as a functions of a transverse field and J 2 
for lattice sizes A'^ = 100 and Ji = 1. 



FIG. 7: (Color online) Three-dimension of the ground state 
fidelity as a functions of magnetic field and Hamiltonian pa- 
rameter J2 for iV = 100 and 5h = 0.001 




d^g/dh 



FIG. 6: (Color online) Three-dimensional plot of the deriva- 
tive of geometric phase versus h and J2 for A'^ = 100 and 
Ji = -1. 




FIG. 8: (Color online) The ground state fidelity versus h and 
J2. The parameters set as A'^ = 100, 5h — 0.001 and Ji — — 1. 



and fidelity in thermodynamics limit specifies that the 
QPTs lines that come from the fideUty approach is con- 
sistent with the geometric phase results which we do not 
show here. 



V. UNIVERSALITY AND SCALING OF 
GEOMETRIC PHASE 

To better understand the properties of the geometric 
phase, and the relation between GP and quantum criti- 
cality, in this section we investigate the scaling behavior 
of GP by the finite size scaling approach. In Fig. ([9| the 
derivative of /3g with respect to the magnetic field has 



been shown for different system sizes for Ji = 1, J2 — 2. 
In thermodynamics Hmit d(3g/dh diverges as the critical 
point is touched, while there is no divergence for finite 
lattice sizes. As the size of system becomes large, the 
derivative of GP tends to diverge close to the critical 
point. More information can be obtained when the max- 
imum values of each plot and their positions are exam- 
ined. 

As it manifests the divergences of d/3g/dh occur at 
— I and h^^ — \/3 where exactly correspond to the 
critical points that were obtained using the energy gap 
analysis {hc^ = /i,r,^c2 = ^o)- The position of the first 
maximum point {JiMaxx) of dPg/dh tends toward the crit- 
ical point like hMax^ = h^, + N'^^ (6*1 = 1.72 ± 0.01) 
which has been plotted in Fig. ( [lO| . Moreover, we 
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FIG. 9: (Color online) The derivative of Pg with respect to 
magnetic field as a function of h for various system size for 
Jl = 1, J2 = 2. 




; h =h 



FIG. 10: (Color online) The scaling behavior of hMaxi in 
terms of system size (A*'), where hMaxi is the position of first 
maximum in Fig.([9|. Inset: Scaling of the position hMax2 oi 
dPg/dh for different-length chains (Ji — 1, J2 ~ 2). 



have shown that the position of the second maximum 
{hMax2) of dj3g/dh goes to the second critical point, such 



as h 



MaX2 



with 6*2 ~ 1.72 ±0.01 (Fig. (10 1, 



inset). 

A more detailed analysis manifest the scaling behavior 
of dPg/dh at the maximum points versus N. In this way 
we have plotted the scaling behavior of df3g/dh\fij^^^^ [i — 



1, 2) in Fig. ( 11 ) and inset of Fig. ( 11 ) respectively, which 
show a linear behavior of df3g/dh\iij^j^^, versus In(A^) with 
the exponents ~ 0.32 ± 0.01. 

To study the scaling behavior of j3g around the critical 
points, we perform finite-scaling analysis, since the max- 
imum value of derivative of l3g scales logarithmic. Then, 



FIG. 11: (Color online) The logarithm of the firs maximum 
of d/Bg/dh versus the logarithm of chain size, In(Ai'), which is 
linear and shows a scaling behavior (Ji — 1,J2 ~ 2). Inset: 
The scaling behavior of the second maximum point of dPg/dh 
for different-length chain in region (II). 




FIG. 12: (Color online) The finite-size scaling analysis for 
the case of logarithmic divergence around the first maximum 
point {hMaxi) for Ji — 1, J2 = 2. The GP, considered as 
a function of system size and coupling, collapses on a sin- 
gle curve for different lattice sizes. Inset: A manifestation 
of finite-size scaling of dPg/dh around the second maximum 
point {JiMax^) for various system sizes in region (II). 



by choosing a proper scaling function and taking into ac- 
count the distance of the maximum of Pg from the critical 
point, it is possible to make all the data for the value of 

*5^P[^~ ^^AJaxi] ^ function of N^/" [h-hMax,) 
for different N collapse onto a single curve. The mani- 
festation of the finite-size scaling is shown in Fig, 
and its inset for several typical lattice sizes. It is c 
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that the different curves which are resemblance of vari- 



7 




h 



N=I2 
N=I6 




FIG. 13: (Color online) The derivative of GP versus h in 
region (I) for Ji = 1, Ja = 0.8. 



FIG. 15: (Color online) Finite-size scaling of d/Bg/dh for dif- 
ferent lattice sizes in region (I). The curves which correspond 
to different system sizes clearly collapse on a single curve. 
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FIG. 14: (Color online) Scaling of the position hmax of d/Sg/dh 
for different-length chains for Ji = 1, J2 ~ 0.8. Inset; Scaling 
of the maximum of d(3g/dh for systems of various sizes in 
region (I). 



[huax) versus In(Af) has been plotted in inset of Fig. 
(14 1. The exponent for this behavior is r ~ 0.32 ± 0.01. 



We illustrate the finite-size scaling behaviors of jSg around 
its maximum points in Fig. (15). It shows that the GP 



can be approximately collapsed to a single curve. These 
results show that all the key ingredients of the finite- 
size scaling are present in these cases too. In these cases 
scaling is fulfilled with the critical exponent = 1, in 
agreement with the previous results and the universality 
hypothesis. We have also investigated the GP behavior 
for Ji = — 1, J2 < 1. Our calculations show that the 
non-analytic and scaling behavior of GP are the same 
as the former results with the same finite size scaling 
{v — 1). The relation between the ground state fidelity 
and Berry phase has been studied in terms of geometric 
tensors [30| . A similar analysis show that the universality 
and scaling behavior of the ground state fidelity are the 
same as geometric phase counterparts in this model. 



VI. SUMMARY AND CONCLUSIONS 



ous system sizes collapse to a single universal curve as 
expected from the finite size scaling ansatz. Our result 
shows that v — 1 is exactly correspond to the correla- 
tion length exponent of Ising model in transverse field 

(^ = 1)- 

However, we have plotted the derivative of fig for Ji — 
1, J2 = 0.8. versus h in Fig. (13 1 for different lattice 
sizes which shows the singular behavior as the size of the 
system becomes large. As it manifests the divergences of 
d[3g occurs at = The similar analysis shows the 

scaling behavior of the position of the maximum point 
(hMax) tends toward the critical point like hMax = he — 
iV"^ (e ~ 1.72 ± 0.01) which has been plotted in Fig. 

dh 



(14 1. The linear behavior of -jS- at the maximum point 



In this work we have studied the geometric phase and 
ground fidelity of the one-dimensional extended quantum 
compass model in a transverse field. We use the Jordan- 
Wigner transformation to construct an explicit analytic 
expressions of geometric phase and ground state of the 
model. We show how the geometry phase and fidelity 
how could detect the quantum phase transitions in the 
inhomogeneous system. Moreover, we have investigated 
the universality and scaling properties of derivatives of 
geometric phase to confirm the results were obtained us- 
ing the correlation functions analysis. The results show 
that the derivatives of the geometric phase diverge close 
to the critical points and exhibit beautiful scaling law. 
Comparing the scaling behavior of geometric phase and 



correlation functions TH' shows that the scaling behavior the following form 
of GP is the same as scaling behavior of correlation func- 
tions, and obtained exponent {v = 1) and universality 
behaviors (logarithmic divergence) confirm our previous 
results. Finally, there is a good agreement between the 
theoretical prediction and numerical data [29] , 
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VII. APPENDIX 



A. Unitary Transformation 



The unitary transformation matrix U which can trans- 
form the Hamiltonian Eq. (fT])into a diagonal form, has where 



J' J* 




)1 




2 _ 2 







\-L-{2h±E^) „ 


((2/1 T E^) (i( J*2 - L*') - L*{2h ± E^f) 


)i 






) - L*{2hi 




~\m-jL*' 





J*{ja-\L\^)- JL*^ 



a — q,p, ^ — 1 for a — q and 7 — 1 for a = p. 



B. Ground State 



By using the unitary transformation the unnornialized 
ground state is obtained as 



\G{h)) = Y[[v,\0) + V2 cfc'^^m + V, 4.VJ,|0) + V, c'J.cf |0) + V, c^JjO) + v, cf c^^^cf c^JjO)], 

k 
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_ 6Ah^ - J^L*^ + J*^(2J^ - + 2| JpC - IQEgh"^ - 4/i(| Jp - |Lp + C)Eg + 8/t^(3| - 2|L|^ + 2C) 

~ 8/l2|J|2+2|J|2C + 2|J|4- J2L*2_ J»2i2 

Eci'ih- Eg){LJ* - JL*)/2 

- 8/l2| J|2 + 2| J|2C + 2| J|4 - J2i*2 _ J*2^2 ' 

_ (4/i, - EG){Ah^J + J( J*^ - i.*'^) + J*C) 

~ 8/l2| J|2 + 2| J|2C7 + 2| J|4 - J2L*2 _ 7*2^2 

_ {Ah - EG){^h^J + J*{.P ~ L^) + JC) 

~ ~8/l2| J|2 + 2| J|2C7 + 2| J|4 - J2L*2 _ j,2_^2 

I 

in which C = \/o? — h. 
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